Geometrical frustration of an extended Hubbard diamond chain in the quasi-atomic 
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We study the geometrical frustration of extended Hubbard model on diamond chain, where vertical 
lines correspond to the hopping and repulsive Coulomb interaction terms between sites, while the 
rest of them represent only the Coulomb repulsion one. The phase diagrams at zero temperature 
show quite curious phases: five types of frustrated states and four types of non-frustrated ones, 
ordered antiferromagnetically. Although decoration transformation was derived to spin coupling 
systems, this approach can be applied to electron coupling ones. Thus the extended Hubbard model 
can be mapped onto another effective extended Hubbard model in atomic limit with additional three 
and four-body couplings. This effective model is solved exactly through transfer matrix method. 
In addition, using the exact solution of this model we discuss several thermodynamic properties 
away from half filled band, such as chemical potential behavior, electronic density, entropy, where 
we study the geometrical frustration, consequently we investigate the specific heat as well. 
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I. INTRODUCTION 

The Hubbard model is one of the simplest model for 
strongly interacting electron systems. In general, rigor- 
ous analysis of the model is a difficult task, only in par- 
ticular case a number of rigorous exact results have been 
obtained^. On the other hand, geometrical frustration 
in strongly correlated electron systems have attracted a 
great deal of interest over the past decades^* 3 -. For in- 
stance, due to the competition between the nearest and 
the next-nearest exchange coupling, the inorganic spin- 
Peierls compound exhibits a transition from a gapless 
phase to a gaped dimerized ground state*. The quan- 
tum phase transition point from the gapless phase to the 
gapped dimerized phase of this model was first deter- 
mined by Okamoto and Nomura 5 . 

The interplay between geometrical frustration and 
strong electron correlation results in a complicated phase 
diagram, containing many interesting phases: spin- 
gaped metallic phase, disorder magnetic insulation phase, 
Heisenberg insulator The Hubbard model on the tri- 
angular lattice shows a transition from a paramagnetic 
metal to an antiferromagnetic insulator as the Hubbard 
on-site repulsion gradually grows 7 . 

Using the numerical diagonalization of finite size sys- 
tem Hida predicted the appearance of the magnetization 
plateau in one of the pioneering theoretical papers^. In 
the framework of the transfer-matrix and dynamical re- 
cursive approaches the frustrated magnetization plateaus 
were obtained for ferromagnetic-ferromagnetic-antiferro- 
magnetic, kagome chains as well as zigzag ladder with 
multi-spin exchanges^—. 

For small nanoclusters by the exact numerical diag- 
onalization the average electron density, magnetization 
plateaus via chemical potential or magnetic field have 
been studied in Hubbard model— 

Recently geometrical frustration of Hubbard model 



was widely studied, particularly the diamond chain struc- 
ture such as considered by Derzhko et aL 14 i 15 , where 
they discuss the frustration for a special class of lattice. 
Montenegro-Filho and Coutinho-Filhoi^ also considered 
the doped AB2 Hubbard chain, both in the weak cou- 
pling and the infinite- U limit (atomic limit). They stud- 
ied a quite interesting phases as a function of hole doping 
away from half-filled band. Mancin i 17 ! 18 has presented 
the exact solution of extended one-dimensional Hubbard 
model in the atomic limit, obtaining the chemical poten- 
tial plateaus of the particle density, of the on-site poten- 
tial at zero temperature and studied the thermodynamic 
charge susceptibility, compressibility at finite tempera- 
ture, as well as other physical quantities. Vidal et alJ^ 
also discussed two interacting particles evolving in a dia- 
mond chain structure embedded in a magnetic field. As 
the particles interact, it leads to the strong localization 
induced by the magnetic field for the particular value 
of a flux. Analogous model was also studied by Rossler 
and Mainemer— . The Hubbard model in other quasi- 
one-dimensional triangular structure also was studied by 
Wang 21 . The latter indicated that for small hopping 
term, the system exhibits short range antiferromagnetic 
correlation, whereas, when the hopping terms become 
greater than critical point, there is a transition from an 
antiferromagnetic state to a frustrated one. Moreover, 
the insulator-metal transition takes place at hopping in- 
teraction even greater than another critical point. Fur- 
ther Gulacsi et aL 22 ' 23 also discussed the diamond Hub- 
bard chain in a magnetic field and a wide range of prop- 
erties such as flat-band ferromagnetism and correlation- 
induced metallic, half-metallic process. The spinless ver- 
sions of the Hubbard model on diamond chain was also 
recently investigated through exact analytical solution 24 , 
as well as Lopes and Dias^ 5 - performed a detailed inves- 
tigation using the exact diagonalization approach. 

In the last decade several diamond chain structures 
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were discussed. Some of them were motivated by real ma- 
terials such as (^3(003)2(011)2 known as azurite, which 
is an interesting quantum antiferromagnetic model, de- 
scribed by Heisenberg model on generalized diamond 
chain. Honecker et ali2& studied the dynamic and ther- 
modynamic properties for this model. Besides, Pereira 
et ali^l investigate the magnetization plateau of delocal- 
ized interstitial spins on diamond chain, as well as they 
detect magnetocaloric effect in kinetically frustrated dia- 
mond chain=L Quite recently, Lisnii 29 studied a distorted 
diamond Ising-Hubbard chain and that model exhibits 
geometrical frustration also. A further investigation re- 
garding the exact evidence for the spontaneous antiferro- 
magnetic long-range order in the two-dimensional hybrid 
model of localized Ising spins and itinerant electrons was 
discussed by Strecka et ah 30 i 31 . Moreover, the thermo- 
dynamics of the Ising- Heisenberg model on diamond-like 
chain was widely discussed in references 3 ^— , and also 
the thermal entanglement was considered by Ananikian 
et alM. 

On the other hand, the analytical exact solution is 
rather amazing, since the exact result is always useful 
to manipulate against the numerical results. Therefore, 
the exact solutions of the models are of great importance, 
so our main goal is to obtain an exact solution for the ex- 
tended Hubbard model on diamond chain in quasi-atomic 
limit. The research of the extended Hubbard diamond 
chain model without the hopping of electrons between the 
nodal sites is based on the three main reasons. First of all 
1/3 magnetization plateau, the double peaks both in the 
magnetic susceptibility and specific heat was observed 
in the experimental measurement o 37 i 38 according to the 
experiments of the natural mineral azurite. Theoreti- 
cal calculations of the one-dimensional Hubbard model, 
as well as the experimental result of the exchange dimer 
(interstitial sites) parameter and their descriptions of the 
various theoretical models are presented. It should be 
noted that the dimers (interstitial sites) exchange much 
more strongly than those nodal sites. There were pro- 
posed various types of theoretical Heisenberg approxi- 
mate methods: the renormalization of the density matrix 
renormalization group of the transfer matrix, the density 
functional theory, the high temperature expansion, Lanc- 
zos diagonalization on a diamond chain are to explain 
the experimental measurements (magnetization plateau 
and the double top) in natural mineral azurite 3 -^. All of 
these theoretical studies are approximate. There is an- 
other possibility. Since dimer interaction is much higher 
than the rest, it can be represented as an exactly solvable 
Ising-Heisenberg model. In addition, experimental data 
on the magnetization plateau coincide with the approx- 
imation of Ising-Heisenberg mode l 27 i 32 ' 36 i 40 and the ex- 
tended Hubbard model without the hopping of electrons 
between the nodal sites 24 on a diamond chain. This is 
the first reason. 

The second one is the quantum block-block entan- 
glement, carried out by exact diagonalization technique 
in one-dimensional extended Hubbard model, calculated 
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Figure 1: (Color online) Schematic representation of 
extended Hubbard model on diamond chain. 



for finite size (L=10) 41 . When the absolute value of 
nearest-neighbor Coulomb interactions (our case) be- 
comes smaller, the effect of the hoping term and the on- 
site interaction cannot be neglected. Finally, we would 
like to point out that although the results, obtained in 
this paper are for the 2-site (dimer) system, their quali- 
tative features are the same as for big size system ones. 

And, the third reason is the experimental observation 
of the double peaks both in the magnetic susceptibil- 
ity and specific hea t 27 i 39 ' 42 may be described exactly by 
the extended Hubbard diamond chain model without the 
hopping of electrons or holes between the nodal sites. 
This phenomenon is particularly important in the quan- 
tum case. 

This paper is organized as follows: in Sec. [H] we present 
the extended Hubbard model on diamond chain. We have 
studied the phase diagram at zero temperature in Sec. IIIII 
Further, in Sec. IIV1 we present the exact solution of the 
model. In Sec. [Vl we have discussed the thermodynamic 
properties of the model, such as electronic density, in- 
ternal energy, compressibility, entropy and specific heat 
away from the half-filled band. Finally, Sec. I VII contains 
the concluding remarks. 



II. THE EXTENDED HUBBARD MODEL 

The extended Hubbard model on diamond chain will 
be considered in this paper, as represented schematically 
in FIG.[TJ In the present model we consider the hopping 
term t between sites a and b, besides, there are onsite 
Coulomb repulsion interaction term denoted by U and 
the nearest neighbor repulsion interaction term, denoted 
by V. We assume also that this model has an arbitrary 
particle density, so, we will include a chemical potential 
term noted by fj,, therefore the Hamiltonian of this model 
can be expressed by: 

N 

H = Y^H ll+1 , (1) 

i=i 

with N being the number of cell (sites a, b and c), whereas 
H iti+1 is given by 
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H iti+ i = -t j2 iL f, '." + a v i, L 

^r— I + ^ ' 



er=J.,t 



-M«+^ + i(n? + n9 +1 )) + 

+ V (ntn\ + i« + nf)(n? + , 

(2) 

with ai )CT , and bi j(T (aj^ and b]^) being the Fermi an- 
nihilation (creation) operator for the Hubbard model, 
while by a we represent the electron spin, and by nf a 
we mean the number operator, with a = {a, 6, c}, using 
this number operator we define also the following opera- 
tors nf = nf t + nf^. 

In order to contract the Hamiltonian ([2]) , we can define 
the following operators 



Pi,i+l ~ 2~( n i + n i+i)> 



(3) 



using these operators, we can rewrite the Hamiltonian 
@, as follows: 



H 



<r=4-,t v 

- (M-Vp iii+1 )(n?+nJ) + 



(4) 



The symmetry particle-hole can be analyzed in a sim- 
ilar way, as discussed in reference^!. The phase diagram 
in the half-filled band, can be obtained using the follow- 
ing restriction for the chemical potential /j, = U/2 + 2V. 
This relation could be obtained by using the symmetry 
particle-hole in a similar way, as discussed in reference^. 



III. THE PHASE DIAGRAM 

In order to study the phase diagram at zero temper- 
ature, first we need to diagonalize the Hamiltonian (0| 
at sites a and b. Where by m = {0, ..,4} we will denote 
the total number of mobile electrons per unit cell at sites 
a and b, while by n c we mean the number of electrons 
per unit cell at site-c. The Hamiltonian at sites a and 
b can be written as 16 x 16 matrix, but this matrix can 
be constructed block matrix, where in the largest 
block matrix we have 4x4, therefore the eigenvalues and 
eigenvectors of this matrix are expressed below as: 



(a) State with m — particle 



An 



Dr 



1 5'<> 



|0,0), 



(5) 



where D m = -{p - Vp iA+1 )m- fip iti+1 + Uq iji+1 , 
While the state \Sq) = |0,0) means there are no 
particle at site a and b, respectively. 



(b) State with m = 1 particle 

\W=Dx±t, |5( ± )) = ^(|0,a) T | ( x,0», (6) 

with |S'ct ± '') we mean there is one particle at site a 
or b with either spin up or down. 

(c) State with m = 2 particles 

\ 2 * = D 2 + V, \S 2 „} = W,a), (7) 
Aj } = D 2 + U, \S^) = ^ (||t, 0) - |0,|t)) , (8) 

X^=D 2 + V + 6 -f, (9) 

\sf ] ) = £ (|#,o) + let) + e T (| u> + 1 l,t») (io) 

\f=D 2 + V, |sj t 2) ) = ^(|t,l>-|l,t)), (11) 
where 



C±=V 2 + ?, (12) 

o ± = u-v±>/(u-v)> + ™? t (13) 

the states \S2a-) an d ) are defined in a similar 
way as defined in the previous case, two particles 
both with spin up or down and two particles with 
opposite spins. 

(d) State with m = 3 particles 

a£ } = D 3 + U + 2V ± t, \S$) = ^ (||t, a) T kit)) , 

(14) 

the states isj^) correspond to two particles with 
opposite spin and one spin 



(e) State with m = 4 particles 



A 



D 4 + 2U + 4V, 1 5 



4-1 14/ 



lit, It)- 



(15) 



the state ISItit) rneans there are two spins at each 
site with opposite spins. 

It is worth to note that, the Hamiltonian ([2]), has 256 
eigenvalues per diamond plaquette. 

From here after, we will consider only repulsive on- 
site Coulomb interaction (U > ) as well as repulsive 
Coulomb interaction between nearest neighbor (V > 0). 
Besides V < U, in order to make more realistic our dis- 
cussion. 

Afterwards we will discuss the phase diagram at zero 
temperature, the present model exhibits 9 states which 
are tabulated in table HI for all possible numbers of par- 
ticles available in the chain. For the empty particle 
and fully filled particles, we have one state \Sq) = |Vac) 
and | Se) = (Full), respectively. While for one particle 
per unit cell and five particles per unit cell, we also 
have one corresponding frustrated state | Si) = | FRUi) 
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m 


Ground state energy 


Ground state 








|S ) =|Vac) = Ui\So}i ® |0)i 


1 


—fj, — t 


|Si> =|frui) =u i \sir ) )i® \0)i 


2 


-2fi + V + %i 


|S 2 ) =|AFM 2 ) = IL |S| t " ; )i® |0>i 


2 


-2fi-t + V 


|S 2 )=|FRU 2 )=n i |5«- ) ) 1 ®|l), 


3 


ft t 

-3fi + 3V + %i 


|S 3 ) =|FRU 3 > = IU l4F } )i ® |1>< 


4 


-4ju + 51/ + L7 + ^ 


|S 4 ) =|AFM 4 ) =n i |5'j t - ) )i® |2)i 


4 


-4/i - 1 + 5V + 1/ 


|S 4 )=|FRU 4 )=n i l^- ) ) l ®|l) s 


5 


-5ju 4- 8V 4- 2(7 - t 


|S 5 )=|FRU 5 )=n i |5^ ) )^l 2 )> 


6 


-6p + 12V + 3U 


|S 6 )=|Full)=n i l%4t>i®|2)i 



Table I: The first column n* represents the number of 
particles per unit cell, in the second and third column is 
tabulated the ground state energy per unit cell and its 
corresponding ground state respectively. 




0.2 0.4 0.6 

v/u 



and | S5) = IFRU5), respectively. However, for 2 parti- 
cles per unit cell, we have two states, one configuration 
is an antiferromagnetic state ( | S2) = IAFM2)), whereas 
other configuration is a frustrated state ( | S2 ) = IFRU2)). 
Due to particle-hole symmetry the analyze for the case 
of four particles (or two holes) becomes analogous for 
the case of two particles, hence, we have one antiferro- 
magnetically ordered state ( | S4) = IAFM4)) and another 
frustrated state ( | S4) = IFRU4)). Furthermore, for the 
special case of half-filled particles (with 3 particles or 
holes per unit cell) we also could have two possible states, 
one frustrated state | S3) = IFRU3) and another antiferro- 
magnetically ordered state ( | S3} = IAFM3)), however the 
case IAFM3) occurs only when V > U, this case we ig- 
nore, since, it is too artificial one. In FIG. Owe illustrate 
the phase transition at zero temperature of these states. 
In FIG.[2fa) we illustrate the phase diagram of V/U ver- 
sus n/U, for fixed value of t/U = 1, where we display 
seven states Si (i = 0..6). While in FIG.FJIb) we display 
the phase diagram of t/U versus n/U , for fixed value of 
V/U = 0.25. This phase diagram illustrates 9 states, 
the seven states are already shown in FIG. FJJa) and two 
additional ones IFRU2) and IFRU4). The fully transi- 
tion boundary functions between states are tabulated in 
table fll] The first and third columns correspond to the 
boundary states, while the second and fourth columns 
correspond to the boundary functions. 

In FIG. by dashed line we represent the half-filled 
band curve of the extended Hubbard model on diamond 
chain. 

Other variants of the phase diagram at zero tempera- 
ture could also be discussed, although the main proper- 
ties already had been discussed in FIG. [2j 



IV. EXACT SOLUTION 

The method to be used will be the decoration trans- 
formation early proposed by Syozi^ and Fisher—, more 



Figure 2: Phase diagrams at zero temperature: (a) 
displays the phase diagram of V/U versus fi/U for fixed 
value ott/U = 1. Whereas in (b) displays the phase 
diagram of t/U versus fi/U for fixed value of 
V/U = 0.25. 



State 


V x n/U (fig. 2a) 


State 


t' x n/U (fig. 2b) 


So, Si 


-1 


So, Si 


-t' 


Si, S2 


3+V'- v / (V r '-l) 2 +16 


Si,S 2 


11— 10-\/0.81+16t' 2 +20t' 


2 




S2, S3 


2V' 


S2, S3 


1 


S3, S 4 


2V' + 1 


S3, S4 


6 


S4, S5 


-l+7V'-^/(V"-l) 2 +16 


S4, S5 


17-20t' +5\/0. 81 + 16t' 2 


2 


2(1 


S5, S6 


W + 2 


S5, S6 


t'+i 






S5, S4 


13 
10 






S4, S3 


-it' +3+2^/0. 81+16i' 2 


4 






s 3 ,s 2 


20t' + 13-10V n -81 + 16t' 2 


9X1 






S2, Si 


1 

in 



Table II: The first and third columns represent the 
boundary between two states tabulated in table U while 

the second and fourth columns mean the function of 
boundaries. For simplicity we use the following notation 
V = V/U and t! = t/U 



later this approach was the subject of study by Ro- 
jas et al4^, for the case of multi-spins. As well as by 
Streckai£ for the hybrid system (e.g. Ising-Heisenberg) . 
Another interesting variant of this approach is also dis- 
cussed by us, where we propose a direct transformation 
instead of several step by step one^i, an illustrative suc- 
cessful application of the last approach was performed in 
reference 48 . The decoration transformation approach is 
widely used to solve spin models, however the decoration 
transformation approach can also be applied to electron 
coupling system, such have been used by the authors 24 



5 



for the case of spinless fermion on diamond structure. 

Therefore, in order to use the decoration transforma- 
tion approach we can write the Boltzmann factor for ex- 
tended Hubbard model on diamond chain as follows: 



w c a - c-P Do 



+e~ 
+e~ 



-pn 



( 



e~ pu + 3e 



where the 9± already were defined in equation (1131) . 

The extended Hubbard model on diamond chain can be 
mapped onto an effective extended Hubbard-like model 
in atomic limit, where it involves additional three-body 
and four-body interaction. Therefore, the effective ex- 
tended Hubbard model becomes like: 



H 



Un h n ti 



n i I "'H J.. 



(17) 



where /2 could be interpreted as effective chemical poten- 
tial, in a similar way, U labels the onsite Coulomb repul- 
sion coupling, meanwhile V corresponds to the nearest 
neighbor repulsion coupling, the next terms, W3 will be 
interpreted as three body interaction term, whereas W4 
as four body coupling. All the above parameters could 
be obtained as a function of Hamiltonian @ parameters, 
when decoration transformation is performed 

The Boltzmann factor for effective Hubbard model in 
atomic limit becomes like: 



where the Boltzmann factor is given by Eq. (fl"B| . and 
using a convenient notation, such as 



w , (x) = 1 + 2.T 1 + ^^ 
V z V 
3 1 1 



— + r 

T 
yz 



(21) 



y 4 z &1 



with x — e 2 ^, y = e^ 1317 , z 



e=^* and 

<^ = e'^ \/(u-v) 2 +i6t 2 ^ Whereas the remaining Boltz- 
mann factors could be expressed easily through the func- 
tion wo, 0(1) defined in Eq. (|21[) as follows: 



x (n 1 +n 2 )/2 
y [n 1 /2} + [n 2 /2} 
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(22) 



by [. . .] inside braces we mean the least integer of any 
real number. Thereafter, all elements of matrix (|2TJ)) are 
well expressed just in terms of wo.o{x). 

The Boltzmann factors for the effective Hubbard model 
(I18p with three and four-body terms are given by: 



w ,o 
u>o,i 

U>0,2 
U>1,1 
Wl,2 
W 2 ,2 



A, 

Ar; r = e^ 2 , 
Ars; s = e"^/ 2 , 
ArH- t = e-&, 
Ar 3 st 2 u; u = e ~^ 2 , 
Ar 4 s 2 t 4 u 4 v. 



(23) 



The determinant of transfer matrix becomes a quartic 
equation of the form det(T- A) = A 4 +a 3 A 3 +a 2 A 2 +a 1 A, 
where the coefficients become: 



Ae~ 



3Hi, 



(18) 



the factor A is an additional variable to be determined in 
terms of the original Hamiltonian ^ . By the use of dec- 
oration transformations^— we will impose the following 
condition 



«i 



a 2 



-2^0,0^1,1^2,2 + 2Wq,2 w 1,1 + 2^0,1^2,2 + 
-2w 0i0 wl 2 - 4^0,2^0,1^1,2, 



~2w 2 A + 2100,0^1,1 + »o,ow 2 ,2- 
-1-2101,1102,2 - W0.2 - 2w lj2l 
-ioo,o - w 2 , 2 - 2l0i,i, 



(24) 



(19) 



using Eq. (flU)) we can find all parameters of the effective 
Hamiltonian (| 1 T[) and the factor A. 

To solve the effective Hubbard model with up to four- 
body coupling, we can use the transfer matrix method 49 , 
similarly to that used in referenc o 24 ' 50 . Therefore we 
symmetrize the Hamiltonian by exchanging i — > i + 1 and 
i + 1 — > i, thus, the transfer matrix becomes symmetric, 
for our case, this transfer matrix could be expressed by: 



WOfi WQ,l W0,1 W ,2 

W ,l 101,1 1Di,l 1Dl,2 

WQ,1 101,1 Wis 1Dl,2 

11)0,2 11)1,2 10l,2 11)2,2 



(20) 



or alternatively the coefficients also can be expressed, 
using the effective Hamiltonian parameters, thus, the co- 
efficients of the quartic equation are given by: 



ai = A (-2r 6 < 5 s 2 uS + 2r 6 s 2 t + 2r e s 2 t 4 u 4 v+ 

+2r 6 s 2 t 4 u 2 - 4r 6 s 2 t 2 u) , 
a 2 = A (~2r 2 + 2r 2 t + r 4 s 2 t 4 u 4 v + 2r e t 5 s 2 u 4 v- 

-r 4 s 2 - 2r 6 s 2 t 4 u 2 ) , 
a 3 = -A (l + r 4 s 2 t 4 u 4 v + 2r 2 t) . 

Therefore the roots of algebraic quartic equation may 
be reduced to a cubic one, solutions of which are given 
in the following way: 

A i = 2 v ^co S (^)-ia 3 j — 0, 1, 2, (26) 
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with 



R = 



90203 — 27ai — 2ag 
54 ' 



(27) 



Furthermore, we also have additional trivial solution 
A3 = 0, of the algebraic quartic equation. 

Hence, the largest eigenvalue of the transfer matrix will 
be Ao, which is expressed by Eq. (f2T>|) . Once known the 
largest eigenvalue of transfer matrix, we are able to ob- 
tain the partition function straightforwardly, and there- 
after the free energy is given by / = — ^ln(A ). Using 
the free energy per unit cell, we may obtain the thermo- 
dynamic properties and how the model behaves when the 
number particles are changing away from the half filled 
band. 



V. THERMODYNAMIC PROPERTIES 

In order to study the thermodynamic properties we 
will use the exact free energy / = — — ln(Ao) as a start- 
ing point. Therefore, we will proceed our discussion of 
thermodynamic properties as a function of temperature, 
chemical potential and electronic density. Assuming that 
we are considering only repulsive onsite Coulomb inter- 
action (U > ) as well as repulsive Coulomb interaction 
between nearest neighbor (V > 0). Besides V < U, in 
order to make more realistic our discussion. 



A. The electronic density 



We will explore the electronic density p 



dp' 



as 



a function of chemical potential as well as the hopping 
term. In FIG. El we plot the chemical potential p/U 
versus t/U, for fixed value of temperature T/U = 0.01 
and nearest neighbor coupling V/U — 0.1. With gray 
scale, we mean electronic density between < p < 2, 
the darkest region corresponds to the lowest electronic 
density, while the brightest region corresponds to fully 
filled electrons. 

In FIG. IU(a), we plot the electronic density as a func- 
tion of chemical potential at low temperature for fixed 
value of hopping term t/U = 0.5 and V/U = 0.1, show- 
ing six plateaus at p = 0, 1/3, 2/3, 1, 4/3, 5/3 and 2, 
this phenomenon vanishes as soon as the temperature in- 
creases. Whereas in FIG.Sfb) we plot the same quantity 
but for large hopping term t/U = 2 and V/U = 0.1. In 
this case the plateaus corresponding to electronic densi- 
ties p — 1/3, 2/3, 4/3 and 5/3, which turn away from 
the half filled band p = 1. Moreover, the plateaus for 




Figure 3: Electronic density per site for T/U = 0.01 and 
V/U = 0.1, as a function of the t/U and p/U, where the 
black region corresponds to empty particles, while the 
white region corresponds to fully filled electrons or 
empty holes, by different levels of gray regions we 
represent the intermediate electronic density. 



densities 2/3 and 4/3 become larger, as we can show in 
FIG. |H which is also in agreement with FIG. [3l 



B. The internal energy 



The internal energy S — 



dln(Ao) 

dp 



, for extended Hub- 



bard model will be discussed as a function of electronic 
density. In FIG. [3 a) the internal energy is plotted for 
small hopping term t/U = 0.2 and V/U = 0.1, this 
internal energy exhibits a gap energy (step like func- 
tion) at zero temperature, when the electronic density 
changes, but when temperature increases the step like 
function becomes smoother. Whereas in FIG. EJb) for 
hopping term t/U — 0.5 and V/U = 0.1 is plotted and 
the shape of curves behaves basically similar plot to that 
t/U = 0.2. In summary, FIG. [5] displays the electronic 
density dependence as well as temperature dependence. 
We can see as far as the larger is the electronic density, 
the lower becomes the energy, although for low density 
around p — 0.5 a maximum of internal energy. We also 
notice for several fixed values of temperature, the plot 
for the internal energy becomes smooth for sufficiently 
temperature. 



C. Compressibility 

Another interesting amount that we could discuss 
will be the compressibility defined by n(p) = 
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Figure 4: (Color online) Electronic density as a function 
of chemical potential, for fixed value of V/U = 0.1. (a) 
For t/U = 0.5 and in (b) for t/U = 2.0 
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Figure 5: (Color online) Internal energy <?(p) as a 
function of electronic density p for given value of 
V/U = 0.1. (a) For i/17 = 0.2 and in (b) for t/U = 0.5 



which is convenient to study electron particles, in- 
stead of that compressibility sometimes known as total 

compressibility^ kt = — Af-^ = — ■ Therefore the 

compressibility will be discussed as a function of Hamilto- 
nian parameters, temperature and electronic density. In 
FIG. [HI we display the compressibility behavior against 
the electronic density on the diamond chain, for fixed 
value of V/U = 0.1. In FIG. HJa) we fixed the hopping 
term at t/U — 0.2, roughly below the phase transition 
occurs at t/U = 0.3, where we display how the com- 
pressibility behaves at low temperature. The compress- 
ibility becomes minimum (harder) close to the following 
fractional electronic densities, 0, 1/3, 2/3, 1, 4/3, 5/3 
and 2. Whereas there are 6 local maximum (softer) be- 
tween fractional densities 1/3, 2/3, 1, 4/3, 5/3 and 2, 
for temperature lower than T/U < 0.05. Whereas for 



0.05 < T/U < 0.5, there is just one minimum at p = 1, 
besides the minimum at p — and p = 2. Furthermore, 
for higher temperature T/U > 0.5 the compressibility 
exhibits just one just one maximum at electronic density 
p = 1. In FIG.IUJb) we illustrate for t/U — 0.5, roughly 
above the phase transition. In principle, the behavior 
is quite similar to FIG. Ela), a low compressibility sur- 
vives for higher temperature t/U = 0.1 at half filled band 
p = 1. Although, when temperature is sufficiently high 
(0.3 < T/U < 0.5) there are only one minimum at p = 1, 
besides empty (fully) electrons, whereas for T/U > 0.5, 
showing just a simple maximum at p = 1. 
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Figure 6: (Color online) The compressibility n(p) as a 
function of electronic density p, for a fixed values of 
V/U = 0.1. (a) For t/U = 0.2, and in (b) for t/U = 0.5. 



D. The entropy 
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Figure 7: Entropy as a function of t/U and n/U, for 
fixed value of U/V = 0.1. Black region corresponds to 

non-frustrated entropy while the white region 
correspond to higher residual entropy effect, gray states 
correspond to intermediate residual entropy effect. 



Furthermore, we will study another interesting quan- 
df 

tity the entropy S = — , and how the entropy behaves 

when the Hamiltonian parameters, temperature or even 
electronic density changes. In FIG. [7] we illustrate the 
magnitude of entropy, as a function of t/U and (J./U, as- 
suming a fixed value of V/U — 0.1. By different levels 
of gray scale, we represent the magnitude of entropy, the 
darkest region corresponds to lowest entropy while the 
brightest region corresponds to the higher entropy. At 
low temperature, the large value of entropy is related to 
the influence of the residual entropy. As far the tem- 
perature increases, the entropy curves become soften, as 
displayed in FIG. Eta) for T/U = 0.01 and in FIG. Efb) 
for T/U = 0.05. 

Mean while in FIG. [HI we display the entropy <S(/i) as a 



function of chemical potential, assuming a fixed value of 
V/U = 0.25, where we are able to illustrate the residual 
entropy effect, in agreement to FIG. [7] 

The entropy indicates several peaks via the chemical 
potential. In FIG. [H[a) we plot for t/U = 0.2 where 
is shown the influence of two types of residual entropy 
around of S = ln(2) and ln(3). Meanwhile, in FIG. fflb) 
we display for t/U — 0.5. In this case we only display 
the influence of the residual entropy of S = ln(2), cer- 
tainly with the temperature increases we have the in- 
crease of entropy as well, however for large the entropy 
is smaller than for small fj, i=s 1. 

Whereas in FIG. EE we display the entropy as a func- 
tion of electronic density assuming a given parameter 
V/U = 0.1, for fixed value of parameter V/U = 0.1. 
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Figure 8: (Color online) Entropy S against p at low 
temperature assuming a fixed value of V/U = 0.25. (a) 
For t/U = 0.2, and in (b) for t/U = 0.5 



Figure 9: (Color online) Entropy against electronic 
density assuming fixed Coulomb coupling V/U = O.f . 
(a) For t/U = 0.2, and in (b) for t/U = 0.5. 



For high temperature T ^ 2, the entropy reaches its 
highest point at p = 1. We can note the residual en- 
tropy is strongly related to the electronic density. We 
also indicate, there is no residual entropy for a fully elec- 
tronic density therefore the entropy becomes lower even 
for high temperature. 

For electronic density p — 1/3 and 2/3 the residual 
entropy is S = In (2), while for p = 4/3 and 5/3 larger 
residual entropy is displayed S = ln(4). 



FIG. \W\ considering a fixed value of V/U = 0.25. In 
the low temperature limit, we can see the effect of the 
phase transition at zero temperature with very sharp 
peaks when the temperature decreases. Particularly, in 
FIG. [TOT al we display the specific heat as a function of 
p/U for t/U = 0.2 while in (b) the specific heat as a 
function of p/U is illustrated for t/U = 0.5. 



E. The specific heat 



d 2 f 

Finally we will discuss the specific heat C = 

behavior for extended Hubbard model in quasi-atomic 
limit. Thus, let us start displaying the specific heat as 
a function of chemical potential, which is illustrated in 



Lastly in FIG. [TT] we display the specific heat as a 
function of temperature assuming the nearest neighbor 
Coulomb interaction V/U = 0.25. In FIG. [TlTa') we con- 
sider for t/U = 0.2, while in FIG. [TlT b) we assume for 
t/U = 0.5. From where we can see very sharp peaks 
as the temperature decreases, this anomalous peaks ap- 
pears due to low lying energy around first order phase 
transition at zero temperature. 
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Figure 10: (Coloronline) The specific heat against Figure 11: (Coloronline) The specific heat versus 

chemical potential for low temperature, for a fixed value temperature, for a given value of V/U = 0.25 and 

of V/U — 0.25. (a) For t/U — 0.2, and in (b) for different values of chemical potential. In (a) is 

t/U = 0.5 considered for t/U = 0.2, while in (b) for t/U = 0.5. 



VI. CONCLUSIONS 

The proposed Hubbard model on diamond chain, was 
discussed at zero temperature as well as at finite tem- 
perature. The phase diagram at zero temperature dis- 
plays 4 frustrated states and 5 non-frustrated states an- 
tiferromagnetically ordered. Concerning to finite tem- 
perature properties, this model can be solved exactly 
through decoration transformatio n 43-45 ' 47 ■ Transfer ma- 
trix technique 4 ^ presents mapping the proposed model 
onto an exactly Hubbard model in atomic limit with 
three and four body couplings. Therefore, detailed ther- 
modynamic properties were discussed, such as density as 
a function of chemical potential, illustrating 6 plateaus. 
Internal energy was also considered as function of elec- 
tronic density far away from the half filled band, showing 
for several fixed temperatures. Besides we studied the 



compressibility. Once again we illustrate this amount as 
a function of electronic density, and so we can conclude 
that the diamond chain is more compressible at low tem- 
perature, when the electronic density is respectively be- 
tween 1/3, 2/3, 1, 4/3, 5/3. As soon as the temperature 
increases, the compressibility decreases, turning into a 
simple decreasing curve as a function of density. We also 
considered the entropy as another interesting amount. 
This amount was studied as a function of chemical po- 
tential as well as electronic density, where we clearly ob- 
serve the residual entropy contribution, owing to geomet- 
ric frustration. Finally we discuss the specific heat, as a 
function of chemical potential and the temperature. In 
the nearest future, we plan to study the Coulomb inter- 
action between the nodal (the monomer-monomer) sites 
and a different type of interaction between the nodal- 
interstitial (dimer-monomer) sites, both positive or neg- 
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ative ones, at the quantum level, which explains the wide 
range of experiments on the diamond chain. 
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